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ABSTRACT 

We present the J- and Kg-band galaxy counts and galaxy colors covering 
750 square arcminutes in the deep AKARI North Ecliptic Pole (NEP) field, us- 
ing the FLoridA Multi-object Imaging Near-ir Grism Observational Spectrometer 
(FLAMINGOS) on the Kitt Peak National Observatory (KPNO) 2.1m telescope. 
The limiting magnitudes with a signal-to-noise ratio of three in the deepest re- 
gions are 21.85 and 20.15 in the J- and X^-bands respectively in the Vega mag- 
nitude system. 

The J- and -K^-band galaxy counts in the AKARI NEP field are broadly in 
good agreement with those of other results in the literature, however we find some 
indication of a change in the galaxy number count slope at J ~ 19.5 and over 
the magnitude range 18.0 < K s < 19.5. We interpret this feature as a change in 
the dominant population at these magnitudes because we also find an associated 
change in the B — K s color distribution at these magnitudes where the number 
of blue samples in the magnitude range 18.5 < K s < 19.5 is significantly larger 
than that of K s < 17.5. 

Subject headings: cosmology: observations — galaxies: evolution — galaxy: formation 
— surveys — infrared: galaxies 
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Introduction 



The infrared universe still holds many unresolved questions: for example, the formation 
and evolution of galaxies and large scale structure. To answer such questions, JAXA 
(Japan Aerospace eX ploration Agency ) launched an infrared astronomy satellite AKARI 



(formerly ASTRO-F, 



Murakami 



20041 ) on 22 February, 2006 (Japan Standard Time, JST). 



AKARI is equipped with a 68.5 cm c ooled telescope and two scientific instruments 



namely the 



Onaka et al- 



ar- In frared Surveyor (FIS, iKawada et al.ll2004l ) and the Infrared Camera (IRC 



2004 ) . A major goal of the mission is to perform an all-sky survey in six 



infrared wavebands from 9-180 microns to higher sensitivity, spatial resolution and larger 
wavele ngth coverage than the fi rst such survey made by the Infrared Astronomical Satellite 



(IRAS, 



Hacking fc Houck 



19871 ). In addition AKARI will perform deep surveys over the 
wavelength range of 1.7 microns (near-infrared) to 180 microns (far- infrared) . Since the 
orbit of AKARI is Sun-synchronous polar, AKARI can often obtain deep exposures at 
the ecliptic poles. Such a deep survey program c urrently underway in the AKARI North 



Ecliptic Pole (NEP) field flMatsuhara et al 



20061 ). Our near- infrared (NIR) ground survey 



is designed to identify the counterparts of sources detected with the AKARI NEP deep 
survey. We will study the nature of the detected sources via examination of their broadband 
energy distributions by also using the optical deep Subaru images in the same field (Wada 
et al. 2007 in preparation). In this paper, we mainly discuss the J- and Kg-band galaxy 
counts and compare these with the results from other fields. 

The galaxy counts as a function of magnitude are a basic quantity in the investigation 
of galaxy evolution. It has been studied extensively in the optical region. However the NIR 
galaxy counts provide more accurate evolutional information on galaxies, since the effects 
of dust extinction and star formation are less severe than in the optical, hence the NIR 
light traces the underlying stellar mass more closely. Furthermore if-corrections due to the 
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redsh ifted spectra 



type (jCowie et al. 



energy distributions are also smaller and nearly independent of galaxy 



1996 : 



Poggianti 



19971 ) 



For over a decade, there have been numerous surveys in the NIR i-T-band. The results 
of these galaxy counts have stimulated a lot of debates. The first claim against the single 
power-low slope of the .fT-band galaxy counts was reported by I Gardner et al.l (119931 ). From 
four blank-field surveys (HWS, HMWS, HMDS, HDS), they showed that the slope of their 
counts changed at K ~ 17 from d(logN)/dm ~ 0.67 to 0.26. In co ntrast, subsequent 
deeper counts with the Keck 10 m telescope ( Djorgovski et al.lll995l ) did not show any 
evidence for a turnover down to a limiting magnitude of K ~ 24. However a superposition 



of the available results for wide and deep counts indeed showed a 



at 16 < K < 20 (jSaracco et al. 



1993) 



Cristobal-Hornillos et al. 



jump in the count slope 



(120031 ) point out that 



reproducing the change in the slope at K$ ~ 17.5 requires a delay in the galaxy formation 
epoch for both elliptical and spiral galaxies to lower values, Zf orm < 2 and the presence 
of a population of star-forming dwarfs at all redshifts in a A-dominated universe. The 
tendency of the bump is also seen in recent result s from wide and de ep surveys, though the 



fainter counts of wide field surveys (e.g. 7.1 deg 2 



Elston et al. 



20061 ) re quire a correction 



for completeness w 



Caputi et al 



2005 



lile the brighter part o 



Forster Schreiber et al 



coun ts of deep surveys (e.g. 



Roche et al. 



2003 



20061 1 have large errors due to the lower number 



.fT-band count slope tend to have bluer optical-N] 


R colors (e.g. 


Gardner et al. 


1993; 


Gardner 


1995: 


Diorexrvski et al. 


1995: 


Saracco et a . 


1997: 


Bershadv et al. 


19 


98; 


Saracco et al. 


1999; 


McCracken et al. 


2000; 


Huang et al. 


2001; 


Kiimmel & Wagner 


2001 


). On the other hand, 



papers describing the J-band galaxy counts 



published to date nu mber only a few. The 



deep counts with the Keck 10 m telesco pe (Bershadv et 
covering a total area of roughly 1 deg 2 (jVaisanen et al. 



al. 



2000 



998) and wide field surveys 



Drory et al 



20011) did not 



show any sign of a flattening in the slope down to their limiting magnitudes. It is important 
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to examine the existence any such change of galaxy counts and corresponding color by 
using the data obtained from improved deep, wide-field imaging surveys, since these will be 
fundamental clues to understand the evolution in the galaxy populations. 

The i^s-band survey in the AKARI NEP field covers the magnitude range 
16.5 < Kg < 19.5, so is therefore useful for determining the exact location of the change in 
galaxy count slope and the change of galaxy color. On the other hand, the J-band survey 
covers the magnitude range 17.0 < J < 20.5, and is therefore valuable to link other deep 
and wide galaxy counts. We also make an attempt to discover any corresponding change or 
confirm the unreported change in the J-band galaxy count slope by using our larger galaxy 
sample. 

We describe the observations and reduction strategy in §|2J We then present our 
procedure to separate stars and galaxies, completeness corrections and completeness 
corrected galaxy counts in §31 The galaxy color trend and the color distribution are describe 
in §H In §5} we discuss the source of the uncertainties in the galaxy counts, possible 
causes of observed change in the galaxy counts and the bluer trend in the B — K$ color 
distribution. Finally, the summary of this research is given in §6j Throughout this paper 
a cosmology with H = 70 kms _1 Mpc _1 , Q m = 0-3 and Q\ = 0.7 is assumed in ag reement 



with the results for Type la supernovae (IRiess et al 



1998 



Perlmutter et al 



WMAP flSpergel et al 



19991 ) and 



20031 ). Unless otherwise noted, we use the Vega magnitude system. 



2. Observations and Data 
2.1. Observations 



The J- and Ks-bsmd imaging were made with the FLAMINGOS (lElston et al.lll998 . 

20031 ) on the KPNO 2.1 m telescope in 2004 June 12-16. The detector is a Hawaii II 
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2048x2048 HgCdTe science grade array divided into four quadrants with eight amplifiers 
each. The FOV is 20.7' with pixel scale of 0.606". To cover the entire field of the optical 
Suprime-CAM images (of which details are given in Wada et al. 2007 in preparation, four 
different pointings, which are named as NE, NW, SE and SW, were used in each band. 
The total effective area is 750 arcmin 2 in both the J- and i^s-bands. Figure [T] show these 
divided survey areas on the z'-band image of the AKARI NEP field. The center coordinate 
and area of these observations are listed in Table [TJ 

The observation in each band and position follows the standard sequence which consists 
of the 5x5 dither pattern with 20"offsets. A single exposure time of 120s was employed 
throughout the observations in the J-band, while the 20s and 30s exposure times in the 
-fTs-band were used to adjust to the background caused by the thermal emission from warm 
optics and telescope structure. The total exposure time in the J-band for each segment 
was 9480s, 8400s, 13080s, and 10920s in the NE, NW, SE and SW, respectively and in the 
Jis-band was 6160s, 7110s, 6950s and 7980s in the NE, NW, SE and SW, respectively. The 
sky condition during this observation was photometric and we also confirmed that the fluxes 
of sources on the data were stable through the observation. The stellar image size of these 
regions ranged from 1.08" to 1.86" in FWHM for the J-band, and fro m 1.08" to 2.04" in 



Persson et al 



(1998) 



FWHM for the i^5-band. The standard stars, 9169 and 9177 listed in 
were observed for flux calibration. Uncertainties in the flux calibration are estimated to be 
less than ±0.1mag. Table [2] summarizes the observational conditions i.e., the total exposure 
time, FWHM of point sources, the 50% completeness limit, and the 3a limiting magnitude 
for each position. 
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2.2. Reduction 



Data reduction was based on the IRAF data reduction package q Dark-averaged 
images for each exposure time were made and subtracted from all images. The sky flat 
images made with a lot of dark sky images were not practical because they were not free 
from any thermal component coming from the warm optics, and telescope structure. Thus 
flat-field calibration images were made by the subtraction of lamp-off dome flat images 
from lamp-illuminated dome flat images. These flat-field images are used to flatten all the 
individual images. 

A extern al IRAF package, t 



(XDIMSUM, 



Stanford et al 



;he Experimental Deep Infrared Mosaicing Software 
19951 ) are used to perform careful sky subtraction, bad-pixel 
and cosmic ray masking, and final combination of the sky-subtracted individual images. We 
describe the details below. 

In the first stage, the 'median' sky image made with five neighboring images were 
used to subtract for the sky component of each image which also consists of the thermal 
component from the warm optics and telescope structure. The possible cosmic ray events 
or rare pixel events due to unstable behavior and tracks of artificial satellites were excluded 
by rejecting the pixel values with scatters exceeding 2a over the entire images. Then, we 
generated the reduced images in each band and each field in the sky, stacking all of the 
sky-subtracted frames with relative offsets which had been determined from the location of 
bright sources. From the first reduced image, we performed source extraction to make a 
object catalog as the input of the second stage reduction. This catalog was used to make 



X IRAF is distributed by the National Optical Astronomy Observatories, which are oper- 
ated by the Association of Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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the object masks for each image. Then we repeated the same steps of the first stage; sky 
subtraction, bad-pixel and cosmic rays masking, and final stacking of the sky-subtracted 
individual images. The difference from the first stage was to make much more accurate sky 
images for the sky subtraction. To make the new 'median' sky image of five neighboring 
images, We masked out the objects some of which could not be recognized on each image. 
Thus we succeeded with an accurate sky subtraction using the new sky images which were 
free from faint objects. Final stacked images of the second stage were the final products in 
this data reduction. 

For source detection and ph otometry on the r e duced images, We employed the 



SExtractor routine, developed by 



Bertin fc Arnoutd (119961 ). Detection was calculated using 



the rms map generated during the image combination process. Objects were detected in the 
condition with four pixels rising more than 1.5a above the sky background rms noise, but 
only those brighter than 3a limiting magnitude in Table [2] were listed in the source catalog. 



3. J- and fCg-Band Galaxy Counts 

3.1. Star-galaxy separation 

Star-galaxy separation is important especially at bright magnitudes, since the modeling 
of the faintest galaxy counts to study the evolution of distant galaxy populations relies 
heavily on an accurate normalization of the models at brighter magnitudes. However the 
spatial resolution of the NIR images is not good enough for a morphological separation of 
stars and galaxies to be secure. Instead, the detected objects were classified on the basis of 
the image profile of their optical counterparts (Wada et al. 2007 in preparation). 

First, we defined visible counterparts within a 2.0" radius from the location of the J- 
or i^5-band selected sources. In our results, the objects with only one optical counterpart 
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for a J- and K$ source are 96.8 and 98.9%, respectively. Objects having two or more 
optical counter parts are rejected from the catalog. We also checked bright counterparts 
reaching to the saturation limit. Figures [2] and [3] show the FWHM of the z'-band image as 
a function of z'-band magnitude. Saturated point-like objects are seen at z' < 18. We found 
that optical counterparts at J ~ 17 and K$ ~ 16.5 reaching to the saturation limit are all 
point-like stars. Therefore we regarded objects below these saturation limits as reliable. 
On the other hand, some faint objects galaxies may be identified as stars because their 
extended emission vanishes into the background noise. However most objects at the faint 
end are visible as extended sources above the background noise (approximately more than 
80% of the total number of the objects at J > 19.5 and K$ > 18.0). The classification into 
point-like and extended objects is therefore reasonably reliable from J ~ 17 and K$ ~ 16.5 
to a 3(7 limiting magnitude listed in Table [2J 

Secondly, we made use of the star-galaxy separation on the color-color diagram to 



obtain an educated guess at the object size. F igure EE shows t 



ae B — %' vs. i' — Kt 



color diagram whic h is similar to the figures of 



McLeod et al. 



Huang et al. 



(119951); 



Oardnerl dl9 95): 



(119971 ) except that we use the i'-band on the SDSS system (IFukugita et al. 



19961 ). B and %' images from SUBARU/Suprime-CAM (Wada et al. 2007 in preparation) 
have been used to derive B — i' and i' — K$ colors using SExtractor AUTO magnitude. The 
limiting magnitudes with a signal-to-noise ratio of five are B ~ 28.4 and %' ~ 27.0 in the 
AB magnitude system, respectively. 

In Fi gure HI the so l id lin e indicates the color of F0-M4 stars calculated based on the 



library of 



Pickles et al. 



(119981 ). The open diamonds represent point-like objects from the 
catalog and flock around the stellar track. The extended objects shown by dots usually 
have a redder i' — Kg color. Therefore, all objects of more than 5.4 pixels (corresponding 



to 1.08") in z' image are regarded to be galaxies. 
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We also checked that our star-galaxy separation was reasonable by comparing other 



star counts in the same region. Figure shows a comparison of point- 
diamonds are our star counts and the open squares represent those of 



ike star counts. The 



Kiimmel & Wa gner 



Kiimmel fc Wagnerl (120001 ) 



(120001 ) . Our star counts at 16.5 < Ks < 17.0 are lower than 
because our star counts are classified on the basis of the optical image profiles and some 
saturated point-like objects are removed from the catalog. However, the fainter star 
counts are in good agreement with e ach other and stellar contamination at K > 18 is less 



important (jKummel &: Wagner 



20001 ) . In fact point-like objects of less than 5.4 pixels are 



approximately less than 20% of the total number in our catalog at Ks > 18 and J > 19.5, 
respectively. 



3.2. Completeness correction 

We calculated completeness correction factors by applying Monte Carlo simulations to 
real images, using the STARLIST and MKOBJECTS packages within IRAF. This approach 
is appropriate because we can take account of the real noise and the pixel-to-pixel sensitivity 
variations in the images since randomly embedded artificial objects within the real images 
are extracted in the same way as that used to generate the real source catalogs. In principle, 
morphological types of the galaxies influence the estimate of completeness for the faint NIR 
galaxies for which galaxy morphologies are unknown. However, due to relatively large seeing 
size almost all the image profiles of the J- and Kg-band sources are much smaller than the 
seeing disk (~ 1", e.g. the effective radii of elliptical bulges are only 0."1 — 0."5 at z ~ 1). 
Therefore, we made artificial objects (mock galaxies) as point-like. Since too many artificial 
objects cause heavy blending among themselves, we generated 25 artificial objects in the 
real image and repeated this procedure 2000 times to obtain a statistically robust number of 
objects for each 0.25 magnitude bin. The radii of the mock galaxies are set to the FWHM 
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of each different pointing field listed in Table [2j The probability of an overlap between real 
and mock galaxies is less than two percent. Objects with fluxes near the detection limit are 
also seriously affected by noise and moreover, sometimes suffered from contamination from 
neighbor objects. Such objects are likely to be affected by flux-boosting and may not be 
extracted with accurate magnitudes. In order to judge the reliability of the source catalog, 
we derived the detection completeness without contamination from flux-boosting. Figures 
|6] and [7J show the ratio of the number of embedded artificial objects in each magnitude bin 
to those that were actually extracted. The NW and SW fields of the AKARI NEP region 
are shallower than the other fields in both the J- and i^s-bands (see Table [2]) and show 
a 20% difference in the completeness levels at J ~ 19.3 and K$ ~ 18.3, respectively (see 
Figures [6] and [7]) . However these areas are less than 30% of the total effective area and 
therefore do not significantly influence the total counts. The detection completeness as a 
function of magnitude are presented in Figures [8] and [91 It is noteworthy that the detection 
completeness is only 50-60% at the faint end, while the inverse of completeness correction 
factors in Figures [6] and [7J are 80-100%. The derived completeness correction factors as well 
as the corrected counts for the J- and i^-bands are listed in Table @] and El respectively. 

It is also important to evaluate the reliability of the completeness correction because 
the correction factor directly affects the faint end counts. For that purpose, we checked 
our completeness correction by comparing them with the completeness corrected shallow 
counts which comprise of a fraction of the co-added frames and th e deep raw counts wh ich 



Vaisanen et al. 



2000). 



comprise of all co-added frames (a similar approach was taken by 
As an example, Figure [10] shows the .fCs-band deep raw counts in the NE+SE regions, 
co-added shallow images at the same region and completeness corrected shallow counts. The 
shallow images comprise of 17% of all co-added frames correspond to K$ ~ 18.5 at the 50% 
completeness limit and are around 0.8mag shallower than those of the all co-added deep 
raw counts. Note that the deep raw counts in the fainter range also require a completeness 
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correction and therefore cannot be directly compared with the corrected shallow counts. 
From Figure [TU], it is found that the faint end of 30% corrected shallow counts at Kg ~ 18.4 
are consistent with the deeper raw counts. 



3.3. J-band galaxy counts 



In this and next subsection, we present the results for the J- and i^^-band galaxy 
counts in the AKARI NEP field. Since the area and depth vary among the four fields 
(see Table |5J), we first obtained the galaxy number counts in each individual field and 
then combined them to obtain the final number counts. Actual numbers of galaxies and 
corrected counts per square degree per magnitude, completeness correction factor and 
effective survey areas in square degrees in the J- and -K^-bands are listed in Table 0] and 
Table O respectively. 



In order to evaluate the slope of the J- band galaxy counts, we fitted power laws of the 



form 



N{mag) =ax i 6 (™*ff-i5) 



(1) 



The best-fit slopes are then calculated in two different domains (J > 19.5 and J < 19.5), 
since the counts in different ranges have different slopes (the cause of the observed change 
in the J-band galaxy counts at ~ 19.5 is discussed in 15. 3p . The derived parameters (a, 
b) of the fits for our and the available J-band galaxy counts are listed in table [5J The 
slope of our ga l axy co unts G?(log N)/dm ~ 0.39 at J > 19.5 is in good agreement with 



Vaisanen et al 



(120001 ) whereas the der ived slope of o ur ga laxy counts d()og N) / dm ~ 0.30 



at J < 19.5 favo r the h igher counts of 



of 



Saracco et al 



Teplitz et al. 



(Il999l ) as opposed to the lower counts 



(jl9991 ). The J-band galaxy counts in the AKARI NEP field show a clear 



change in the slope at J ~ 19.5 from 0.39 to 0.30. 
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We also show no-evolution models in both the J- and i^s-bands for comparison. 
The no-evolution models are pure i^-correction models and constructed from three 
galaxy types (early-type, late- type 1 and late-type 2). The ^-corrections of each galaxy 
ty pe are calculated by using the theoretical spectral energy distribution (SED) models 



of 



Bruzual fc Chariot 



(120031 ) . Regarding the if^-ba nd galaxy counts , we r efer to the 



type-dependent K-band luminosity function (LF) of iKochanek et al 



2MASS Second Incr e ment al Release Catalog of Extended Sources (IJarrett et al 



however 



Cole et al 



(2001) based on the 



(1200 ll ) point out that the K-hand photometry of 



2000) 



Loveday et al. 



20001 ) has better signal-to-noise and resolution than the 2MASS images and enables more 



accu rate 2MASS mag nitudes to be measured than the original default magnitudes (see 



also 



Frith et al. 



20051 ). We therefore normalized the characteristic ma gnitude Mt of the 



no-evolution model at K s < 13.25 of the 2MA SS counts measur ed by 



Frith et al. 



Cole et al. 



(120051 ) 



(120011 ). On the other hand, to 



using a similar magnitude estimator to that of 
derive the J-band galaxy count s, we assume the shape of the type-dependent J-band LF is 



the same as the i^-band LF of 



M* K to Mj based on the 



Kochanek et 



Bruzual fe Chariot 



al 



(200ll) and simply convert the normalized 



(120031 ) SEDs. The basic characteristics of the 



SEDs i.e., galaxy type, metalicity and star formation rate are listed with the parameters of 
the LF of each galaxy type in Table [3l The redshift of galaxy formation used for all galaxy 
types is z form = 4. 

For the J-band galaxy counts, the no-evol ution model is consistent with the results 
of the Deep Near-Infrared Sky Sur yey (DENIS 



the Deep Multicolor Survey (DMS 



Martini 



VI anion 



19981 ) and the bright counts of 



20011 ) at magnitudes J < 17. However the 



no-evolution model is far below our observed counts from J ~ 17 to fainter magnitudes; 
the factor is 2.1 at J ~ 19.6. This is because neither pure luminosity evolution (PLE) nor 
bulk (luminosity or number density) evolution in the extragalactic population is included 
in the no-evolution model. In particular, the inclusion of PLE into the no-evolution model 
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will increase the normalization of the models in the intermediate magnitude range but we 
reserve this work (and the investigation of bulk evolution in the near-infrared population as 
a whole) to the next paper in this series. 



3.4. fCs-band galaxy counts 



The .fTs-band galaxy counts are presented in Figure [12l which shows the results in the 
AKARI NEP field along with a compilation of published galaxy counts. In the bright 



magni t ude range Ks < 17.0 . our re sult shows s 



(119931 ): iKiimmel k Wagneri (1200(1 ): 



Elston et al. 



i ghtly 



(I2006h : 



higher counts than 



Gardner et al. 



6.7. The counts are however consistent with those of 



toll); 



Vaisanen et al 



(120001 ) 



Cristobal- Hornillos et al 



by a factor of < 1.5 at Ks 



McLeod et al 



(119951 ): 



Martini 



(120031 ). On the other hand, 



though our faintest counts are approximately 50% complete and therefore subject to large 
correcti ons, they are in goo d agreement with the results of other .fT-band deep surveys 



(except 



Saracco et al. 



19991 1 . The K s -band galaxy counts in the AKARI NEP field show 
an indication of a possible 'bump' in the magnitude range 18.0 < Ks < 19.5. We discuss 
uncertainties and a possible cause of the observed change of the slope in the J- and i^^-band 
counts in §5J 

For the i^-band galaxy counts, although the no-evolution model is consistent with the 
counts derived from the wide surveys for Ks < 17, they are lower (approximately 35%) than 
the present work at around Ks ~ 18. In addition, the no-evolution model cannot reproduce 
the characteristic change in the slope since the count slope of each galaxy component in the 
no-evolution model varies rather smoothly from the bright to the faint end. 
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4. B - K s Galaxy Colors 



4.1. B — Ks galaxy color trend 



Figure [13] is a plot of the B — Ks color trend of the galaxies in the AKARI NEP 
field. The filled circles show the median color in each magnitude bin. All of the galaxies are 
detected at the 3<r limiting magnitude in the -fTs-band listed in Table [2J The median B — Ks 
color of the galaxies becomes slightly redder up to Ks ~ 17.5 where it reaches a maximum 
value of B — Ks ~ 6, and then tends to be bluer to fainter magn i tudes. This trend is similar 



to previous results (j Gardner et al. 



1993 



McCracken et al. 



2000 



Huang et al.l 



200 lh . 



4.2. B — Ks galaxy color distributions 

Figure [TH shows the B — Ks galaxy color distributions. The solid line is the bright 
galaxy sample selected in 16.5 < Kg < 17.5, dashed and dotted lines correspond to the 
galaxy samples selected in 17.5 < K s < 18.5 and 18.5 < K s < 19.5, respectively. These 
samples cover the bump region of the Ks-bsmd galaxy counts described in §3.41 The 
number of sources in the bright sample (16.5 < Ks < 17.5) increases from B — Ks ~ 4 to 
a peak at B — Ks ~ 6.5 and decreases thereafter. The decrease of the red end is steeper 
than the increase at the blue end. Though the peak of middle sample (17.5 < Ks < 18.5) 
is somewhat shifted toward a bluer B — Ks color of 6.2 and the number of blue and 
red galaxies are relatively larger than the bright sample, the shape of the B — Ks color 
distribution of the middle sample is similar to that of the bright sample. On the contrary, 
the shape of the faint sample (18.5 < Ks < 19.5) is obviously different from the others, 
with a dramatic increase in the number of blue galaxies at 3.5 < B — Ks < 5.0. We confirm 
this difference in the two samples by applying a Kolmogorov-Smirnov test (K-S test) to 
our data. The K-S test results in a probability of 10~ 18 that two (bright and faint) galaxy 
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samples are identical. 



5. Discussion 

We presented the change of the J- and i^^-band galaxy count slopes and the bluer 
trend in the B — K$ color distribution in §3.31 §3.41 and §U respectively. In this section, 
we discuss the relation ship of J- and .fCs-band count slopes ( 15. ip . the source of the 
uncertainties in the galaxy counts (15. 2p and implications from other observational studies 



5.1. Relation of J- and i^-band galaxy counts slopes 



In the J-band, we presented galaxy counts based on the largest sample to date which 
bridges the previous deep and shallow galaxy survey counts. We showed the clear change of 
the slope at J ~ 19.5 from d{\ogN)/dm = 0.39 to 0.30. However the change in the slope 
of the J-band counts is less prominent than that in the i^s-band counts where it is visibly 
seen as a 'bump'. 

The brighter component of the K^-b and counts (above the 'bump') is dominated by 



early-type galaxies obeying the PLE (e.g. 



Metcalfe et al. 



1991 



1996 



Pozzetti et al. 



1996J). 



The less prominent slope change in the J-band counts is due to the 



early-type ga 



Saracco et al. 



axies 



(1999) 



having large J — K$ colors (i.e. fainter J -band magnituc e s). In deed, 
reported J — K > 1.5 at K s < 18, while iBershady et all (119981 ) discussed that early-type 
population show J — K > 2 at z > 1. A J — K color of ~1.5 would imply to a corresponding 
turnover in the J-band counts at J ~19.5. 
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5.2. Uncertainties in the Ks-hand count 



For the i^s-band galaxy counts in the AKARI NEP field, we have shown an indication 
of a change in the slope over the magnitude range 18 < K$ < 19.5. Here we discuss 
the uncertainties in the -fTg-band counts over three different magnitude ranges : 18 > Ks, 



18 < Ks < 19, 19 < Kg, and compare in particular with the FLAMEX counts (jElston et al. 



20061 ). a similar NIR survey covering 7.1 deg 2 using the same instrument (FLAMINGOS). 



In the bright magnitude range (Ks < 18.0), the slope of fitted power law of equation CD 
is d(logN)/dm = 0.32 ± 0.06. Our result shows higher counts than the results of FLAMEX. 
The difference between the two surveys is a factor 1.3 at Ks ~ 16.7. This may be attributed 
to the fact that the two surveys adopted a different approach for the star-galaxy separation. 
Star-galaxy separation of FLAMEX was performed by using color-color diagrams while in 
the case of this work, the separation is based on the image profiles of optical counterparts. 
Therefore some stars with profiles of more than 5.4 pixels may be included in the galaxy 
counts. However star counts in th e NEP field do not heavily affect the galaxy counts in the 



fainter magnitude range, K > 18 (jKummel &: Wagner 



20001 ). In fact point-like objects of 



less than 5.4 pixels are less than 20% of the total number in our catalog at Ks > 18. Thus 
the observed turn over from Ks ~ 17 is considered an intrinsic feature. 

In the intermediate magnitude range (18.0 < Ks < 19.0), the slope of our galaxy 
counts is d(log N)/dm = 0.11 ± 0.02. Our f^s-band counts show a flattening and are in 



good agreement with the FLAMEX results. The ten dency of the 



the results of FLAMEX and other deep surveys (e.g. 



Roche et al. 



iattening is also seen in 



2003 



Forster Schreiber et al 



20061 ) 



Caputi et al 



2005 



In the faint magnitude range (Ks > 19.0), the slope of our galaxy counts is 
d(\ogN)/dm = 0.32 ± 0.06. Our faintest counts may be significantly affected by 
incompleteness (approximately 50% complete) and show lower counts than the results of 
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FLAMEX; the factor is 1.3 at Ks ~ 19.3 



other if-band deep survey s 



Cristobal-Hornillos et al. 



loweyer, they are consistent with the results of 



Gardner et al. 



1993 



McLeod et al. 



1995 



Totani et al 



2001 



20031 ). We therefore cannot rule out the possibility that some of 
the reason for the change of slope observed over the magnitude range 18.0 < Ks < 19.5 
may be due to incompleteness since the completeness corrections are non-negligible at these 



magnitudes. 



5.3. Evidence for number density evolution in the early-type population 



To reproduce the turn over in the slope at Kg ~ 17.5, ICristobal-Hornillos et all (120031 ) 
delayed the galaxy formation epoch for elliptical and spiral galaxies to quite low redshift, 
Zform < 2 and demanded the presence of a population of star- forming dwarfs at all redshifts 
in a A-dominated univer se. However, since numerous early-t ype galaxies at z > 2 have 



recently been found (e.g. 



Daddi et al. 



2005 



Labbe et al 



20051 ). it is unrealistic to completely 



extinguish early- type population at z > 2. 

Our result of the B — Ks color distribution also shows the number of blue galaxies 
dramatically increases below the 'bump' to fainter Kg magnitudes. These results indicate 
that the contribution of early-type galaxies to the counts is relatively low in the fainter 
magnitude range, a res ult that has been suggested by several studies including galaxy 



morphologies at z > 1. 



Kajisawa fc Yamadal (120011 ) reported that the comoving density of 



My < —20 galaxies decreases at z > 1 in each morphological category based on the result 
of the HST WFPC2 /NICMOS archival data of the HDF-N. In particular, they showed 
that the number density of the bulge-dominated galaxies conspicuously decreases to 1/5 
between the 0.5 < z < 1.0 and 1.0 < z < 1.5 bins. Other morphol ogical studies also show 



the lack of massiye elli pticals at redshifts greater than z ~ 1.5 (e.g. 



Franceschini et al. 



1998 



Rodighiero et al 



200 lh . 
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Moreover, a similar decrease in the number of red galaxies at higher redshift is indicated 
from studies of extremely red objects (EROs). The EROs are selected to have very red 
optical-to-near infrared colors such as R — K > 5 — 7, I — K > A — 6. The number counts of 
EROs using a sample (K < 21) of 158 EROs in the ELAIS N2 field shows lower counts than 
those predicted by models in which all E/SO galaxies evolve according to PLE. However, 
the counts can be fit by a model combining PLE with galaxy mergin g and a decrease i n 



2002 1. 



the comoving number density of these passive galaxies with redshift (IRoche et al. 
Based o n a very deep samp le (K s < 22) of 198 EROs from the ESO/GOODS data in the 



CDF-S. 



Roche et al. 



(120031 ) reported that the ERO number count flattens at K > 19.5 from 



d(\ogN)/dm = 0.59 to 0.16. 

This negative density evolution also has be en reported from studie s of the evoluti on 



of the .K'-band LF at various redshi 



20031 ). z < 2J 



Bolzonella et al 



and z < 4 (jSaracco et al 



2002 



ts: z < 1 ( 



G 



Pozzetti et al 



20061 ). Note that 



azebrook et al. 



Pozzetti et al 



995: 



Feulner et al 



2003), z < 3. 5 (IKashikawa et al 



20031 ) 



(120031 ) show a contrary result, 



however, their large errors due to the limited sample size do not rule out a decrease of the 
comoving number density with redshift. 

To summarize, the number density evolution of the elliptical galaxy population remains 
unsolved. However the bump in the iT^-band galaxy count is an important clue and imposes 
tight constraints on the formation and evolution of each galaxy population because galaxy 
counts are measured with less indefiniteness than other observables and wider and deeper 
NIR survey can pin down the slope of the galaxy counts more precisely. We entrust further 
discussion about the evolution of each galaxy population to our forthcoming paper. 
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6. Summary 

We performed a NIR survey in the AKARI NEP field, using FLAMINGOS on the 
KPNO 2.1m telescope. The total area is 750 arcmin 2 and the limiting magnitudes to a 
signal-to-noise ratio of three in the deepest regions are 21.85 and 20.15 in the J- and 
.fCs-bands respectively in the Vega magnitude system. The extracted objects are classified 
into stars and galaxies by using image profiles of their optical counterparts from Subaru 
Suprime-cam observations. The resultant star counts are consistent with other published 
results in the same region. Completeness corrections are estimated by adding artificial 
sources to the images and extracting them. The reliability of the completeness corrections 
are also evaluated by comparing the completeness corrected shallow counts and the deep 
raw counts. We then presented the J- and i^-band galaxy counts and the B — Ks color 
distribution in three different magnitude slices down to K s = 19.5. Our J- and K s -b&nd 
galaxy counts are broadly in agreement with those of other published results, while our 
counts indicate a change in the slope in both the J- and -K^-band counts. The slope of our 
J-band galaxy counts changes at J ~ 19.5 from d(\ogN)/dm = 0.39 ± 0.02 to 0.30 ± 0.03, 
while the i^^-band galaxy counts show a more significant 'bump' in the magnitude range 
18.0 < K s < 19.5. The B — K s color distribution of our galaxy samples also shows an 
obvious change around the 'bump' magnitude of the X^-band counts. The number of blue 
galaxies in the magnitude range 18.5 < K s < 19.5 is significantly larger than that of the 
Kg < 17.5 samples. 

The authors wish to thank the referee who's comments greatly improved the contents 
of this work. 
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Fig. 1. — Final coverage for our NIR survey in the AKARI NEP field. The four segments 
marked by solid square boxes with white lines are overlaid on the z'-band image (Wada et 
al. 2007 in preparation). The top-left, top-right, bottom-left and bottom-right square boxes 
are referred to as NE, NW, SE and SW, respectively. The center coordinate and area of 
these observations are listed in Table [TJ 
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Fig. 2. — z'-band FWHM vs. magnitude with the corresponding J-band objects. The 
point-like objects at the seeing limit cluster together around a FWHM of 5.2 pixels. 




Fig. 3. — z'-band FWHM vs. magnitude with the corresponding Ks-baiad objects. Same as 
Figure [2] but for the iT^-band. 
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Fig. 4. — B — i' vs. z' — .fTs color diagram. The open diamonds and dots represent point-like 
objects of FWHM less than 5.4 pixels and extended objects of FWHM more than 5.4 pixels, 
respectively. The solid line shows the main sequence of F0-M4 stars. 
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Fig. 5. — Comparison of star counts in the NEP field. The filled diamonds and the open 



Kummel &: Wagnerl (120001 ) . respectively. 



squares represent the star counts of this work and 
Error bars are calculated from the number of stars in each magnitude bin using Poisson V N 
counting statistics. 
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Fig. 6. — The completeness correction factor, derived from Monte Carlo simulations (see 
text in §3.2! for details) are given for each magnitude bin. Error bars are estimated from the 
standard deviation of the set of the Monte Carlo simulations. 
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Fig. 7. — The completeness correction factor for the i^g-band. Symbols are the same as 
Figure El 



-30- 




17 18 19 20 21 22 

K s [mag] 



Fig. 8. — The J-band detection completeness estimated from Monte Carlo simulations. Each 
line shows the fraction of detected sources in the J-band for each pointing region. 
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Fig. 9. — The i^-band detection completeness. Meaning of the lines are the same as that 
in Figure El 
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Fig. 10. — Comparison of number counts to two different depths. The open diamonds 
indicate the raw counts derived from all co-added frames in the NE and SE regions. The 
filled squares and the open circles show raw counts and completeness corrected counts of a 
subset of the NE+SE data, for which only 17% of all frames are co-added. 



-32 - 



10 5 



D) 
(0 

E 

73 



1 1 r 

Maihara '01 SDF i— ©— 
Martini '01 DMS i— b— 
Saracco '99 NTT >- -o — 
Teplitz '99 : 




15 16 17 



18 
/ [mag] 



Total 

Early-type 

Late-type 1 ■ 

Late-type 2 

This work NEP ♦ 

Vaisanen '00 ELAIS N1 s a 

Vaisanen '00 ELAIS N2 i v- 

Mamon '98 DENIS i — ±- 
i i i 

19 20 21 



22 



Fig. 11. — J-band galaxy counts. The black diamonds are galaxy counts derived from 
this work. Error bars are calculated from the number of galaxies in each magnitude bin 
using Poisson y/~N counting statistics. Note that at the faintest magnitudes, completeness 
errors may also have a non-negligi ble effect on the sour ce counts. The symb ols in the figur e 



correspond to, Maihara '01 SDF ( Maihara et al. 



Vaisanen '00 ELA 



19991), Teplitz '99 (ITeplitz et al. 



S Nl and N 



2 ( Vaisanen et al. 



2001), Martini '01 DMS (Martini 



2000), Sarac co '99 N TT 



1999) and Mamon '98 DENIS (IMamon 



2001). 



Saracco et al. 



19981 ). The counts 



having large errors (more than 50% in relative errors) are not shown. The symbols and 
colors for each survey are the same for the Ks-band galaxy counts. Orange, green and blue 
lines indicate the no-evolution predictions for each galaxy type, described in §3.31 Red line 
corresponds to the total galaxy counts of these predictions. 
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Fig. 12. — fTg-band galaxy counts. The black diamonds are the galaxy counts derived from 
this work. Error bars are calculated from the number of galaxies in each magnitude bin 
using Poisson y/N counting statistics. Note that at the faintest magnitudes, completeness 
errors may also have a non-negligible effect on the source counts. Red and magenta color 
symbols listed in the bottom-right corner are publis hed counts of wi de field surveys. The 



symbols in t he figure correspond to , Frith '05 2MASS (IFrith et al 



Nl and N2 (Vaisanen et al. 



'99 (IMinezaki et al 



2000), Ku mmel '00 NEP ( 



1999), Szoko 



Gardner '96 (I Gardner et al 



y '98 flSzokolv et al 



200 51), Vaisanen '00 ELAIS 



Kummel &: Wagner 



19981) 



19961 ). and Gardner '93 HMWS (IGardner et al. 



2001T). Minezak i 



Huang '97 (|Huang et al 



19970 . 



1993|). Blue and 



cyan color symbols listed in the top-left corner are from pub lished counts of dee p surveys. Th e 



symbols correspond to Totani '01 SDF (Totani et al. 



Saracco '99 NTT (ISaracco et al. 



2003), Minezaki '99 (IMinezaki et al 



19971 ). Moustakas '97 



'93 HDS and HMDS ([Gardner et al 



200ll ). Martini '01 DMS (Martini 



1999), Cristobal-Hornillos '03 ( Cristobal- Hornillos et al. 



2001), 



Moustakas et al 



1999), Saracco '97 E SO KS1 and 



KS2 



19971 ). McLeod '95 (IMcLeod et al 



19931 ). The FLAMEX counts of 



Saracco et al. 



1995) and Gardner 



Elston et al. 



(120061 ) are 



shown as the bold dotted line with the upper and lower errors shown as the corresponding 
upper and lower faint dotted lines, respectively. The counts having large errors (more than 
50% in relative errors) are not shown. The symbols and colors for each survey are the 
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Fig. 13. — B — Ks galaxy color trend for the AKARI NEP survey field. Objects detected in 
both the B- and ii'g-bands are plotted by dots. The filled circles show the median color in 
each magnitude bin. The dash line represents the region of incompleteness; galaxies above 
this line are too faint to be detected in the 5-band (B > 28.5mag in AB). 
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Fig. 14. — B — K s galaxy color distributions for the AKARI NEP survey field. All objects 
are detected in both the B- and i^-bands. Each line corresponds to selected objects in a 
different magnitude bin. 
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Table 1: Observed center coordinates and areas. 



Region name 


R.A. 


Dec. 


Area 




(J2000) 


(J2000) 


(arcmin 2 ) 


(1) 


(2) 


(3) 


(4) 


NE 


17h 55m 48.2s 


+66d 45' 41" 


234 


NW 


17h 53m 41.3s 


+66d 45' 40" 


91 


SE 


17h 55m 47.8s 


+66d 30' 41" 


306 


sw 


17h 53m 42.3s 


+66d 30' 43" 


119 



Note. — (2) Unit of R.A. is hours, minutes, and seconds. (3) Unit of Dec. is degrees, arcminutes, and 
arcseconds. 
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Table 2: Observational conditions of NIR survey in the AKARI NEP field. 



J-band 



Field name 


Exposure time 


FWHM 50% detect 


ion completeness 


3<r limit 




(s) 


(arcsec) (magnitude) 


(magnitude) 


(1) 


(2) 


(3) 


(4) 


(5) 


NE 


9480 


1.38 


20.4 


21.6 


NW 


13080 


1.86 


20.0 


21.2 


SE 


8400 


1.08 


20.5 


21.9 


SW 


10920 


1.50 


20.3 


21.5 


.Kg-band 


NE 


6160 


1.08 


19.2 


19.9 


NW 


6950 


1.80 


18.8 


19.4 


SE 


7110 


1.08 


19.4 


20.2 


SW 


7980 


2.04 


18.8 


19.5 



Note. — (2) Total exposure time for each pointing field. (4) Measurement of completeness is described in 
gSHJ (5) Aperture radius is set a threes of half of FWHM. 
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Table 3: Description of the SEDs, luminosity function parameters and rest-frame colors used 
in the models. 



Galaxy type 


Metalicity 


Star formation rate 


0* 


a 


(J ~ K s ) 




(Z/Z ) 


(Gyr) 


(Mpc^ 3 x 10~ 3 ) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Early-type 


2.5 


Exp. r = 0.5 


1.3 


-0.8 


1.0 


Late-type 1 


1.0 


Exp. t = 1.0 


2.0 


-0.9 


0.9 


Late-type 2 


0.4 


Exp. r = 2.0 


2.0 


-0.9 


0.9 



Note. — (3) The e-folding time of the star formation rate. (4) The parameter <j>* is a normalization 
coefficient that has the dimensions of the number density of galaxies. 
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Table 4: J-band 


raw and corrected galaxy counts. 








J mag 




Completeness 


N 

1 y corr 


Error 


Area 


(magnitude) 


(mag deg z ) 


(mag Meg 2 ) 


(mag-Meg^ 2 ) 


(arcmin 2 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


17.250 


989 


0.98 


1000 


100 


750 


17.750 


1728 


0.96 


1800 


130 


750 


18.250 


2822 


0.95 


2970 


170 


750 


18.750 


4454 


0.94 


4700 


200 


750 


19.250 


5961 


0.91 


6600 


300 


750 


19.625 


7238 


0.87 


8300 


500 


750 


19.875 


7910 


0.84 


9400 


500 


659 


20.125 


9046 


0.75 


12100 


600 


659 


20.375 


8988 


0.66 


13600 


1100 


306 



Note. — (1) J-band magnitude at each magnitude bin. (2) Raw galaxy counts detected in this magnitude 
range. (4) Completeness correction factors. (3) Completeness corrected counts. (5) Errors are dominated by 
Poisson y/~N statistical errors. Errors in the completeness correction factor itself are found to be negligible. 
(6) Total effective area. 
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Table 5: Comparison of the power-law slopes of the J-band galaxy counts. 

J < 19.5 a b Range Area 

(magnitude) (arcdeg 2 ) 

(1) (2) (3) (4) (5) 

NIR AKARI NEP survey 150±30 0.39±0.02 17.0-19.5 2.1 x KT 1 

Vaisanen et al. f2000 . ELAIS Nl) 140±30 0.40±0.03 17.0-19.5 2.0-0.6 x 10" 1 

Vaisanen et al. Q2000, ELAIS N2) 130±30 0.38±0.02 17.0-19.5 5.0-2.9 x 10" 1 
J > 19.5 

NIR AKARI NEP survey 330±140 0.30±0.03 19.5-20.5 2.1-0.9 x 10" 1 

Teplitz et al. (1999) 430±150 0.28±0.03 19.5-21.0 5.3 x 10~ 2 

Saracco et al. (1999 . ESQ NTT) 80±120 0.37±0.11 19.5-21.0 5.5 x 10~ 3 



Note. — (2), (3) The ampliltude of the fitting power law given by eq. [T] (4) Magnitude range of the fitted 
data. (5) Effective area of each survey. 
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Table 6: Kg-hand 


raw and corrected galaxy counts. 








K s mag 


N 


Completeness 


N 

1 11 corr 


Error 


Area 


(magnitude) (mag Meg 2 ) 


(ma 


g -1 deg -2 ) 


(mag _1 deg~ 2 ) 


(arcmin 2 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


16.750 


2707 


0.98 


2760 


160 


750 


17.250 


4406 


0.97 


4500 


200 


750 


17.750 


5702 


0.96 


5900 


200 


750 


18.125 


8025 


0.96 


8400 


400 


750 


18.375 


8908 


0.94 


9400 


400 


750 


18.625 


8390 


0.87 


9600 


500 


750 


18.875 


9067 


0.89 


10200 


600 


540 


19.125 


8427 


0.69 


12200 


700 


540 


19.375 


7435 


0.50 


15000 


1200 


306 



Note. - 



Description of the counts is given in the footnote of Table @] 



